Assembly-oriented reliability analysis method for the top- 


connection structure of a nuclear fuel assembly 


Yingdong Liu*, Wenqiang Li*^*, Jiahao Chen*?^, Sitong Ling?*^, Changfu Wan? 


* School of Mechanical Engineering, Sichuan University, Chengdu 610065, China 
^ Innovation Method and Creative Design Key Laboratory of Sichuan Province, 


Chengdu 610065, China 


Abstract: The nuclear fuel assembly is the core component of a nuclear reactor. In a 
pressurized water reactor fuel assembly, the top-connection structure connects the top 
nozzle to the guide thimble. Its performance reliability 18 essential for the stability of 
the nuclear fuel assembly. In this study, an assembly-oriented reliability analysis 
method for top-connection structures is presented by establishing an assembly-oriented 
top-connection structure parameter modeling method and a nonlinear contact gap and 
penetration correction method. A reliability model of the top-connection assembly 
structure, including multiple stochastic design variables, was constructed, and the 
overall reliability of the top-connection assembly structure was obtained via a Kriging 
model and Monte Carlo simulation. The acquired experimental data were consistent 
with real-world failure conditions, which verified the practicability and feasibility of 


the reliability analysis method proposed in this study. 
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1. Introduction 


Nuclear energy is an important clean energy source [1-3]. The nuclear fuel 


assembly, as the core system in nuclear reactors, governs the safety and reliability of an 
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entire nuclear energy system [4-6]. The top-connection structure is an assembly 
composed of multiple nuclear fuel assembly parts that serve as the connection between 
the top nozzle and the guide thimble in the pressurized water reactor (PWR) nuclear 
fuel assembly. During the lifting and transportation of a fuel assembly, the top- 
connection structure undergoes an acceleration that is a factor of 4 greater than the 
gravitational acceleration [7]. In the case of systematic structural failure, the entire 
nuclear fuel assembly faces significant potential safety hazards. 

To date, reliability analysis of the top-connection structure of nuclear fuel 
assemblies relies heavily on physical experiments, which can be both expensive and 
time-consuming. Numerical simulations have gained popularity in recent years [8-10]. 
However, compared to the numerical simulation reliability analysis of an individual 
part, the analysis of the top-connection structure is a typical reliability analysis of an 
assembly structure that contains several related and mating assembly components. The 
marginal analysis error of any component has a considerable impact on the overall 
assembly reliability. In addition, owing to the existence of different assembly surface 
mating types in the assembly, the gap and penetration relationship of each contact 
surface in the assembly significantly affects the accuracy of the assembly reliability 
analysis. Therefore, developing a reliable analysis method for assemblies remains a 
challenging international focus. 

In recent years, numerous studies have been conducted using numerical 
simulations of the top-connection structure of nuclear fuel assemblies and their 
associated structures [11,12]. Wang et al. [13] established a numerical simulation model 
for a top-connection structure using finite element analysis and investigated the tensile 
process of the top-connection structure. Liu [14] analyzed the strength of the top nozzle 
during transportation and lifting by simplifying the top-connection structure into a fixed 
connection. Yoo et al. [15] established a finite element model of a spacer grid that 
simulated its static and dynamic mechanical properties and compared it with physical 
experiments to verify the accuracy of the model, laying a good foundation for 
subsequent spacer grid optimization. Duan and Zhao [16] utilized CFX software to 


build a geometric model of a bottom nozzle and performed numerical simulations of 


the bottom nozzle at different flow rates and split ratios. Wu et al. [17] conducted a 
mesh sensitivity analysis on the top nozzle, upper core plate, and top grid in the top fuel 
section of a reactor using computational fluid dynamics (CFD and compared the 
calculated pressure distribution with the corresponding experimental results. Su et al. 
[18] investigated the stress and bearing capacity of a CF3 fuel assembly bottom nozzle 
using a finite element analysis and load tests. Their results showed that the stress under 
each condition was in accordance with the American Society of Mechanical Engineers 
(А5МЕ) code and that the bearing capacity satisfied the design requirements of the CF3 
fuel assembly. Wei et al. [19] conducted a numerical study on the coolant flow in the 
top nozzle of a fuel assembly based on CFD and developed a numerical simulation 
method for the top-nozzle drag characteristics. The resistance characteristics were 
analyzed and evaluated to determine the resistance coefficient of the top nozzle. Dyk 
and Zeman [20] proposed a vibration modeling method for guiding the thimble in 
nuclear fuel assemblies and analyzed the effect of the maximum dynamic lateral 
deformation of the guide thimble on the radial gap of the spacer grid. Xu et al. [21] 
developed a CFD model for transverse flow in the top section of the AP1000 core, 
including the top-connection structure, and obtained a reasonable transverse flow 
velocity in the top section of the AP1000 core. Zhao et al. [22] established a nonlinear 
finite element model of a fuel assembly and conducted lateral stiffness, forced vibration, 
and impact experiments to verify its accuracy. 

Despite the aforementioned studies probing the characteristics of top connections 
and associated structures, two major issues remain unaddressed. First, most of the 
studies were simulations of the top-connection structure and associated individual parts, 
while the system reliability analysis of the top-connection structure as an assembly was 
not included. Second, the design parameters in the numerical simulation analysis of the 
top-connection structure were commonly set to constant values. However, under real- 
world working conditions, multiple design parameters, including material properties 
and physical dimensions, have exhibited certain random errors. Hence, the key solution 
is to improve the reliability of the analysis in the presence of random multidesign 


parameter errors. In light of this, in the current study, we propose an assembly-oriented 


top-connection structure reliability analysis method and process with an assembly- 
oriented parametric modeling and correction method for each contact surface in the 
assembly and an assembly-oriented reliability analysis method for the top-connection 
structure. An assembly-oriented approximation model for the reliability of the top- 
connection structure is also established. A Monte Carlo simulation (MCS) was deployed 
to perform a reliability analysis of the top-connection structure with multidesign 
parameter randomness, enabling the failure probability and reliability of the top- 
connection structures to be obtained. The proposed method demonstrates sound 
conformity with the physical experimental data and can be widely adopted in the 


reliability analysis of complex structures with assembly relations. 
2. Assembly modeling and contact surface correction methods 


Unlike the individual parts, an assembly structure is composed of multiple parts. 
Reliability analysis of the assembly system is highly nonlinear because the parts are 
intercoupled, making it difficult to guarantee calculation accuracy. There are three 
crucial elements in safeguarding the accuracy of the reliability analysis of assembly 
structures: determining the parametric modeling of the assembly structure, establishing 
the correlations between each assembly unit, and determining the connection and 
penetration relations between each assembly unit. In this study, a parametric modeling 
method for assemblies and a correction method for the contact gap and penetration of 
the contact surfaces between each assembly part were established, thereby providing 
an effective solution to fluctuations in the accuracy of the assembly reliability analysis. 
2.1 Parametric modeling method for assembly systems 

Assembly parameters are interrelated and mutually restricted. Consequently, the 
correlation between each parameter should be considered to establish the mapping 
relations between every design parameter, as well as the mapping relations between 
design parameters and mating-surface parameters, to ensure that the mating relations 
between the parts remain valid when the assembly model is updated. Mating between 
the parts of the assembly is achieved by matching the features of the parts, and a contact 


surface can be considered as a match of multiple features between two parts. To avoid 


conflicting design variables caused by part modeling when updating the assembly, a 
parametric modeling method for assembly oriented systems was devised. 

First, the design variables of each part were determined, and the corresponding 
features of each design variable were determined accordingly. Subsequently, the 
decoupling relations between the design features and design variables were completed 
by transforming the design-associated matrix to achieve mutual independence of the 
design variables. Therefore, if a certain design variable is changed, it does not cause 
conflicts with other design variables. Finally, the assembly was transformed based on 
features from the contact relations between parts to the contact relations between design 
variables. In this manner, the mating between parts of the assembly can be controlled 


using design variables, as shown in Figure 1. 
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Figure 1. Parametric modeling method for assembly-oriented systems 


The nuclear fuel assembly top-connection structure, as a typical assembly model, 
contains three contact faces: the contact between the top nozzle and connection-sleeve 
outer surface, the contact between the connection-sleeve inner surface and the locking- 
thimble outer surface, and the contact between the guide-thimble and connection-sleeve 
inner surfaces. Parameterization was performed for these contact faces and the design 


variables of each part were determined, as shown in Figure 2. 
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Figure 2. Schematics of the top-connection structure design parameters 

Each part of the top-connection structure contained one parameter, and the design 
variables were independent of each other. Therefore, the mating surfaces of the top- 
connection structure correspond to the design parameters of each part, specifically D1— 
D2, R2-R6, R3-R5, 03-04, апа D6-D3. The links and equation constraints among the 
design parameters were established, as listed in Table 1, to ensure that the mating 
relationships between the parts remained valid when the top-connection structure model 


was updated. 


Table 1. Mapping relations between the design parameters 


R1=D1/2 R2=R3 D2-DI R5-R3 
Н1-Н3-Н2 D4-D3 D6-D3 R4-D3/2 
D3«D2 D5<D4 D7<D6 


Seven part design parameters were obtained by mapping the relations: D1, H2, R2, 
D5, D3, R6, and D7. The remaining eight parameters were obtained from the seven 
parameters by mapping relationships. 
2.2 Contact surface gap and penetration value correction method for assembly 
systems 

Errors occur between the mating surfaces during the assembly process; these have 
direct effects on the gap and penetration values between the contact surfaces. In addition, 


the quality of the modeling mesh affects the contact surface gap and penetration values 
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of an assembly. In this study, a correction method for the initial penetration and gap of 
the contact surfaces was developed. This may be utilized in the following circumstances: 
When there 18 an initial gap between the contact surfaces of the assembly, the gap can 
be supplemented to facilitate surface contact or, in the case of slight penetration 
between the two contact surfaces, the common penetration area can be removed by 
adjusting the contact surfaces. The adjustment method is shown in Figure 3, in which 


the penetration area is marked in red. 


Contact surface gap Gap correction Contact surface penetration Penetration correction 
Figure 3. Contact surface gap and penetration correction method 


The top-connection structure of a nuclear fuel assembly contains three contact 
surfaces, of which two are nonlinear: the contact surface between the top nozzle and 
the connection-sleeve outer surface and that between the connection-sleeve inner 
surface and the outer surface of the locking-thimble outer surface. The other is linear: 
the contact surface between the guide thimble and the inner surface of the connection 
sleeve. Therefore, for the nonlinear contact surfaces of the top-connection structure, the 
gap and penetration values under different mesh sizes can be established by analyzing 
the penetration and gap variation patterns between each contact surface under different 
mesh sizes, as listed in Table 2. 


Table 2. Nonlinear contact surface gap and penetration values for different mesh sizes 


Connection sleeve and locking thimble contact surface Top nozzle and connection sleeve contact surface 


5mm 4mm 3mm 2mm Imm 5mm 4mm 3mm 2mm Imm 
Penetration 
0.45 0.36 0.28 0.0006 0.0001 0.05 0.038 0.034 0.03 0.003 
/mm 
Gap 
1x105 3x10 2x107 1x108 1x10° 3x10? 2х10% 3х10% 2x107 3x10! 
/mm 


As presented in Table 2, the direct penetration value of the contact surface 
decreases as the mesh size between each contact surface gradually decreases. The gap 


between each contact surface also decreases with a decrease in the mesh size. 


Corrections for the gap and penetration values ofthe contact surfaces for different mesh 


sizes were performed, and the results are listed in Table 3. 


Table 3. Penetration and gap correction for each contact surface with different mesh sizes 


Connection sleeve and locking thimble contact Top nozzle and connection sleeve contact 
surface surface 
5mm 4mm 3mm 2mm 1mm 5mm 4mm 3mm 2mm Imm 
Penetration/mm 0.45 0.36 0.28 0.0006 0.0001 0.05 0.038 0.034 0.03 0.003 
Penetration 
2x10 3x10 — 2d0? 9х10"* 3x10 — 2x10" — 2x10? — 8105 — 5x10 — 3x1055 
correction /mm 
Gap/mm 1x10? 3x106 2x107 15109 1x10? 3x106  2x10$ 3х10% 2х107 3x10" 
Gap correction 
0 0 0 0 0 0 0 0 0 0 


/mm 

As presented in Table 3, all corrected gap values were 0. It is assumed that the 
contact surfaces are tightly matched with no gaps. In nonlinear contact analysis, this 
facilitates the complete transmission of force on each contact surface, is beneficial for 
the convergence of the results, and ensures the accuracy of the results. The corrected 
penetration values for different mesh sizes were ай «10 P. Consequently, the 
penetration values between the contact surfaces can be considered negligible. This will 
assist in the adoption of a suitable mesh size to reduce the computational resources 
required for the assembly, as well as converge the analysis to ensure the accuracy of the 


results. 


3. Approximation model construction method for the top-connection 


structure 


Building a reliability approximation model is a prerequisite and foundation for 
reliability analysis of a system. In this study, a method and process for the construction 
of an assembly-oriented reliability approximation model was established using the 
following stages: pre-experimental design of the connection structure, sensitivity 
analysis of the design variables, and verification of the accuracy of the approximation 
model. 

3.1 Pre-experimental design of top-connection structure variables 


In the experimental design for assembly, unlike generic parts, the design variable 


spans of parts tend to overlap. Therefore, conducting a pre-experimental design can not 
only verify the plausibility of multiple design variable spans of the assembly but also 
reduce the number of design variables and preserve those with substantial influence on 
the reliability response. In the mechanical analysis of the top-connection structure of 
the nuclear fuel assembly under transportation and lifting conditions, the loads included 
the weights of the fuel assembly and control-rod assembly and axial acceleration. The 
total weight of the payload was ~7500 N and the axial acceleration was 4 g. Based on 
the ASME code, the third strength theory was used to evaluate the results of this study 
[23]. The simulation results for the top-connection structure under the transport and 


lifting conditions are shown in Figure 4. 
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Figure 4. Stress distribution diagrams of each part 


Figure 4 shows that the overall stress of the assembly under the transportation and 
lifting conditions is mainly distributed near the contact surface. The maximum stresses 
of each part are at the contact surface, being 142.78 MPa for the connection sleeve, 
94.992 MPa for the locking thimble, and 88.78 MPa for the top nozzle. The maximum 
stress value of each part was considered as the output response, and the maximum 
output response of each part was parameterized to realize a closed loop between the 


input and output parameters. The number of top-connection structure design variables 


was finally set to seven: D1, H2, R2, D5, D3, R6, and D7. The output response was set 
as the maximum stress value of each part from the simulation analysis, namely, Stress1, 
Stress2, and Stress3, corresponding to the maximum stress values of the top nozzle, 
locking thimble, and connection sleeve, respectively. 

To ensure a uniform distribution of sample points in the design space and improve 
design efficiency, the Latin hypercube sampling method was deployed for the pre- 
experimental design of the design variables [24]. The design matrix contained 36 
sample points, and the results of the pre-experimental design are presented in Appendix 
1: 

3.2 Sensitivity analysis of top-connection structure design variables 

The number of assembly design variables is typically large, with each design 
variable exerting various degrees of influence on the output response. Therefore, 
reserving all design variables for analysis imposes a considerable burden on 
computational resources. Through a sensitivity analysis, the trend of the influence of 
each design variable on the output response, as well as the specific value of the 
influence, was obtained. To improve the computational efficiency, a sensitivity analysis 
was performed on the assembly design variables, and the design variables with a greater 
influence on the output response were selected to establish a reliability approximation 
model. 

A sensitivity analysis was performed on the design variables of the top-connection 
structure. The extent of influence of design variables D1, H2, R2, D5, D3, R6, and D7 
on the output responses Ѕігеѕѕ1, Stress2, and Stress3 are shown in Figure 5, and the 


specific sensitivity values of each design parameter are listed in Table 4. 
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Figure 5. Extent of the influence of each parameter on the output response 


Figure 5 shows that D1 demonstrated a negative correlation for both Stress1 and 
Stress3 and a positive correlation for Stress2. Among the design parameters, D1, D5, 
and D3 had the greatest influence on Stress3, Stress2, and Stressl, respectively. The 
extent of influence of each design parameter on Stress] was ranked as D3 > D7 > D5 > 
H2 > R2 > R6» D1, the extent of influence on Stress2 was ranked as 05 > D3 > H2 > 


D7 > R6 > D1 > R2, and the extent of influence on Stress3 was ranked as D1 > D5 > 


R6 > H2 > R2 > D7 > D3. 
Table 4. Sensitivity values for each design variable 
: 1 Sensitivity 

ПО Stress] Stress2 Stress3 
DI -22.50 16.35 -100 
H2 -50.43 -57.53 -30.56 
R2 44.04 6.82 18.25 
D5 58.89 100 51.51 
D3 -100 -81.57 -7.13 
R6 41.95 45.71 47.45 
D7 -65.62 -48.13 -16.56 


The construction of the response surface of each parameter versus the output 
response provides a straightforward understanding of the sensitivity of each parameter 
to the output response. Figure 6 shows the response surface of each parameter for the 
output response. From the response surface of each parameter to the output, the 


parameters did not exhibit monotonically increasing or decreasing effects on the output 
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response, and a mutual influence existed between the two design variables. By taking 
H2 and D7 Stress1 as a reference, when D7 = 9 mm, the output response decreases 
monotonically with the increase in H2; when D7 = 10 mm, the output response 
decreases first with the increase in H2 and then increases; and, when D7 = 10.5 mm, 
the output response monotonically increases again with the increase in H2. When H2 = 
4.5 mm, the output response increases monotonically with an increase in D7; when H2 
= 5 mm, the output response decreases and then increases with an increase in D7; when 
H2 = 5.5 mm, the output response increases monotonically with an increase in D7. 
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Figure 6. Response surface of each parameter to the output response 
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Combining the response surfaces of each design parameter for Stress1, Stress2, 
and Stress3 reveals that all seven design variables of the top-connection structure have 
a significant impact on the output response. Therefore, there is no need to scale down 
the design variables. 

3.3 Approximation model construction for top-connection structure 

The pre-experimental design and design variable sensitivity analysis identified the 
design variables that influence the output response. Subsequently, an experiment was 
conducted to identify the design variables and obtain a more accurate approximation 


model of reliability. Appendix 2 presents the experimental design results for the top- 


connection structural variables. 

The 120 generated sample points were divided into two groups, one serving as 
learning points to establish an approximation model and the other as verification points 
to assess the accuracy of the model. The principle of this division was to provide a 
sufficient number of learning points and reserve a certain number of verification points. 
Therefore, the sample points were divided in a ratio of 3:1 [25], with sample points 1— 
90 as the learning points and sample points 91—120 as the verification points. 

The Kriging model, as an efficient interpolation method, can yield not only the 
predicted value of a prediction point but also the variance of that point. Because of this 
advantage, this method was used to construct a reliability approximation model in this 
study. 

Assume that the actual relationship between the response and design variables can 


be expressed as 


W(x) = FG) B*z0), (1) 
where F(x) is a linear combination of polynomial functions, В is a linear regression 


coefficient, and z(x) is a Gaussian stochastic process with a mean of 0. The covariance 


between the two points x; and x; is 
2 А 
Cov[z(x;), 2(х)]=0?А(х,х;0), 0) 
where © is process variance, R(x,,x,;@) is the correlation function at two points x; 


and ху, usually a Gaussian correlation function: 
M 
R(x,,x,,0) =| | exp[-6, G7 —7 Y 1. (3) 
т-1 


where Ө 15 a parameter vector, m is the mth-dimensional vector of the input vector, 
and M is the total dimensionality of the input vector. 0 can be estimated by using the 


maximum likelihood method as 
0 = arg max(-N, In(o?^) - In|R. (4) 
A Kriging model was constructed based on n training points x; (i= 1, 2, ..., п). Let 


Y represent the vector of responses at n training points. The vectors of the regression 


coefficients and process variance can be estimated as 


B=(FTR'FYFTRYY, (5) 


o =F pY R Y-F P), (6) 


r(x) 2 [R(x, x3), R(x,x,;0),..., R(x,x,;0)] , u = F(x)" m (x) UY -F f). (7) 


3.4 Approximation model accuracy verification 

The accuracy of the approximation model is crucial to the accuracy of the 
reliability analysis results, and it can be verified via validation points. The metrics for 
the accuracy verification of the approximation model include the relative maximum 
absolute error (RMAE), root mean square error (RMSE), and coefficient of 


determination R? [26,27]. RMAE is given by 


RMAE = max {| y, - 5, Ї . 


(8) 


where N is the number of verification points, y, isthe true value ofthe ith verification 


point, y, is the predicted value of the ith verification point, and y isthe mean value 


of the verification points. The RMAE was used to characterize the absolute maximum 
residual value relative to the standard deviation of the output value of the sample points. 
The closer the value of RMAE is to 0, the higher the accuracy of the approximation 
model. 


RMSE is given by 


RMSE = (9) 


The RMSE was used to characterize the dispersion of the sample points. The closer the 
RMSE value is to 0, the higher the accuracy of the approximation model. 


Finally, R? is given by 


N 
-3,У 
i=l 


R? Нил are. 
У(Уг-3) 
i=l 


(10) 


R^ was used to characterize the degree of agreement between the predicted and real 


values. The closer the R^ value is to 1, the higher the accuracy of the approximation 
model. 


Figure 7 shows the relationship between the predicted and true values, and the 


errors for each method are listed in Table 5. 
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Figure 7. Kriging goodness-of-fit graphs 
Table 5. Kriging model accuracy validation 
Stress1 Stress2 Stress3 
RMAE 2.1Х 1076 1,2Х 1075 5.8X 10:16 
RMSE 11531079 27751079 1.6X 10716 
R? 1 1 1 


As presented in Table 5, the results verified by each method meet the requirements 
and, consequently, prove that the Kriging model constructed according to this sample 


point meets the accuracy requirements. 
4. Reliability analysis method for the top-connection structure 


The top-connection structure of the nuclear fuel assembly is complex. In practical 
engineering, the stochasticity brought by design variables needs to be considered, as it 
affects the overall structure reliability. The reliability analysis flow for the top- 
connection structure used in this study is shown in Figure 8. The probability distribution 


and the type of probability distribution for each random variable in the top-connection 


structure were determined. Subsequently, the limit-state function of the top-connection 
structure was constructed, and its reliability was calculated via MCS invoking the 
approximation model. This not only addressed the low calculation efficiency of the 


assembly system but also ensured the accuracy of the reliability analysis results. 
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Figure 8. Assembly-oriented reliability analysis flowchart of the top-connection structure 
4.1 Multivariate reliability analysis based on MCS 


i=i+1 


MCS is a robust stochastic simulation method for structural reliability analysis 
[28-30]. It is an excellent contemporary tool for structural reliability analysis because 
its calculation accuracy and convergence speed are not affected by the system 
complexity [31-33]. However, the computational efficiency of this method is limited 
and the required number of sample points is large. In light of this, to address the limited 
computational efficiency of MCS combining the assembly-oriented reliability 
approximation model, a Monte Carlo reliability analysis method based on an 
approximation model is proposed in this study. The specific steps are as follows: 

Step 1: Determine the random variables and designate their probability distribution 
and probability distribution types. 

Step 2: Define the maximum number of MCS simulation runs, N. 


Step 3: Generate a sequence of uniformly distributed random numbers. 


Step 4: Convert the generated random number sequence into corresponding 
random variable values. 

Step 5: Invoke the assembly approximation model to compute the response to the 
current value. 

Step 6: Repeat Steps 3—5 until the maximum number of simulations, N, is reached. 
4.2 Probability distribution statistics of random variables for the top-connection 
structure 

In machining the top-connection structure of nuclear fuel assemblies, the working 
accuracy varies because of the accuracy of the machine tools and instruments, as well 
as the proficiency of the machining operators. Therefore, reliability analysis of the top- 
connection structure requires considering each critical dimension as a random variable. 
Furthermore, the influence of material properties and the unevenness of loading on the 
structural function should also be considered, with the modulus of elasticity E of the 
material and loading force F as random variables. Under normal conditions, the 
dimensionality caused by machining errors conforms to a normal distribution. The part 
dimensions in the top-connection structure as an assembly have upper and lower limits, 
thus conforming to a truncated normal distribution. The modulus of elasticity E of the 
material follows a normal distribution, whereas the force load F fluctuates within a 
certain range, making it conform to a truncated normal distribution. The statistical 
properties of each random variable for the top-connection structure are listed in Table 


6. 


Table 6. Statistical properties of each random variable 
Random Mean  Lowerlimit Upper limit ^ Standard 


ў a Distribution Type 

Variable value value value deviation 

Di/mm 13.45 12.95 13.95 0.289 Truncated normal distribution 
H2/mm 5 4.5 5.5 0.286 Truncated normal distribution 
R2/mm 1 0.75 1.25 0.058 Truncated normal distribution 
D5/mm 11.45 10.95 11.95 0.289 Truncated normal distribution 
D3/mm 12.45 11.95 12.95 0.288 Truncated normal distribution 
R6/mm 0.5 0.45 0.55 0.029 Truncated normal distribution 
D7/mm 11.45 10.95 11.95 0.288 Truncated normal distribution 

E/GPa 199 сэн == 1.1X 104 normal distribution 


F/N 1250 1100 1400 86.243 Truncated normal distribution 


4.3 Limit-state function determination for the top-connection structure 

The structural output response of the top connection must satisfy the ASME code, 
with a membrane stress of «138 MPa and a membrane plus bending stress of «207 MPa. 
The failure mode of the top-connection structure is a typical strength failure, and the 


following structural limit-state function can be established [34]: 


Ё(їх Xp) = sod Fon seed: Ja (11) 
where the variables x,,..,x, are random variables affecting the function, R(x,,..., x, ) 


is the rigidity random variable, and $(x;,...,x,) is the stress random variable. When G > 


0, the structure is in a safe state, whereas, when G < 0, it is in a failure state. 

The top-connection structure, which is an assembly composed of several parts, 
fails when one part fails. Therefore, they can be regarded as a system. For the top- 
connection structure, the specific structure limit-state function is 

G(X) = 207 — max {Stress 1, Stress 2, Stress3} = 0. (12) 

The probability of an individual part failure can be expressed as 


P, = P(GG, x, ) < 0); the probability of the entire top-connection structure failure 
is Р,-5(1-411Б,, where i is the number of system parts and the reliability is 


Б. =1- В. 


4.4 Reliability assessment method for the top-connection structure 

In this study, the MCS method was deployed, invoking the Kriging model to 
calculate the reliability of the top-connection structure. The maximum number of runs 
of the MCS was set to 10,000. The distribution types and data for each random 
parameter are listed in Table 8. The results for each response distribution are shown in 


Figure 9. 
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Figure 9. Cumulative distribution of each response 


Figure 9 shows that the stress distributions of Stress] are all within 130 MPa, 
concentrated in the range of 65-110 MPa, and do not exceed 207 MPa; the stress 
distributions of Stress2 are concentrated in the range of 40-110 MPa, with very few 
exceeding 200 MPa; and the stress distributions of Stress 3 are concentrated between 
100 and 170 MPa, with a few exceeding 200 MPa. The failure probability of each part 


is 
N 
P.=—, 
Дэг (13) 


where N is the number of sample points with G(x) < 0 in the MCS and No is the total 


number of MCS runs. The failure probability of each part can be obtained from Figure 


9 and is listed in Table 7. 
Table 7. Probability of failure and reliability of each part 
Top nozzle Locking thimble ^ Connection sleeve 
Failure probability 0 0.15% 1.85% 
Reliability 100% 99.85% 98.15% 


According to Table 7, the reliabilities of the top nozzle, locking thimble, and 
connection sleeve are 100%, 99.85%, and 98.15%, respectively. This corresponds to a 
real-world failure scenario. The connection sleeve, placed between the top nozzle and 
locking thimble, was directly loaded by the entire fuel assembly during operation. 
Compared to the top nozzle and locking thimble, it exhibited the highest probability of 
failure during operation in the fuel assembly [13]. In general, the structural reliability 
of the top connection remained sound even when considering the randomness of each 


variable. 


5. Conclusion 


An assembly-oriented reliability analysis method for top-connection structures of 
nuclear fuel assemblies was presented. Based on the assembly features, a parametric 
modeling method for assembly systems and a design parameter correlation method 
were proposed, and an assembly-oriented correction method for the nonlinear contact 
surface gap and penetration of the top-connection structure was established to address 
the reliability analysis accuracy fluctuation caused by gaps and penetrations between 
contact surfaces during the assembly analysis. An assembly-oriented reliability analysis 
model for the top-connection structure was established, considering the randomization 
of design variables, such as machining error, material, and force load. The Kriging 
model and MCS method were deployed to analyze the structural reliability of the top- 
connection structure of the nuclear fuel assembly, and the obtained reliability of the 
assembly system demonstrated better agreement than the physical experimental data. 
The proposed method, with sound applicability and practicability, is applicable not only 
to the reliability analysis of the top-connection assembly structure of a nuclear fuel 
assembly but also for the reliability analysis of complex structures with multiple 
assembly relations. Notably, the material characteristics and reliability of nuclear fuel 
assemblies are significantly affected by the long-term operation of nuclear fuel 
assemblies under high temperature, high pressure, and high levels of radiation. Time- 


varying reliability analysis of nuclear fuel assemblies will be the focus of future studies. 


Appendix 1: Pre-experimental design matrix and experiment results 


Table 8. Pre-experimental design matrix and experiment results 


N DI H2 R2 DS D3 R6 D7 Stress 1 Stress2 Stress3 
1 13.263 5.347 1.058 9.688 11.750 0.540 9.896 72.371 51.022 108.638 
2 13.712 5.292 0.914 9.604 11.917 0.476 10.354 63.620 42.801 102.795 
3 12.142 4.875 1.008 9.021 10.917 0.493 9.396 71.481 61.728 119.617 
4 13.487 5.486 0.969 9.979 10.817 0.460 9.271 93.277 124.675 137.641 
5 12.292 4.847 0.964 9.396 11.783 0.549 9.313 68.693 55.227 280.260 
6 14.608 4.736 0.975 9.563 11.450 0.504 9.854 76.737 74.076 124.963 
7 13.413 5.014 0.997 10.479 11.017 0.471 9.021 115.021 222.975 192.132 
8 14.533 4.597 1.075 10.146 11.117 0.518 9.229 101.634 144.597 153.750 
9 12.665 5.181 1.053 9.313 11.083 0.538 9.938 78.716 62.022 109.951 


10 
11 
12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 


12.516 
12.890 
12.441 
13.188 
13.861 
12.815 
14.309 
14.160 
13.637 
12.367 
12.217 
13.114 
14.384 
14.683 
12.591 
13.936 
13.562 
14.085 
13.786 
14.758 
12.964 
14.010 
14.235 
13.338 
14.459 
13.039 
12.740 


5.042 
5.236 
4.986 
4.625 
4.514 
4.569 
5.153 
5.264 
5.069 
4.708 
4.792 
5.208 
5.125 
4.764 
5.458 
4.542 
4.931 
4.903 
5.375 
4.681 
5.319 
4.958 
5.403 
4.819 
5.431 
5.097 
4.653 


1.042 
0.992 
1.031 
1.069 
0.936 
1.036 
0.981 
0.908 
1.086 
0.953 
1.064 
0.919 
0.958 
0.986 
1.014 
0.942 
1.019 
1.047 
1.003 
1.097 
0.947 
0.925 
1.092 
0.931 
1.025 
1.081 
0.903 


10.438 
10.188 
10.229 
9.438 
9.938 
10.104 
9.646 
9.271 
9.813 
10.063 
9.104 
9.188 
10.021 
9.146 
9.354 
9.771 
10.271 
9.896 
9.229 
9.479 
10.396 
9.729 
10.354 
9.521 
9.854 
9.063 
10.313 


11.950 
11.550 
10.950 
11.217 
11.683 
11.717 
11.383 
10.983 
10.850 
11.317 
11.417 
11.583 
11.250 
10.883 
11.483 
11.050 
11.850 
11.883 
11.817 
11.650 
11.150 
11.983 
11.283 
11.183 
11.517 
11.617 
11.350 


0.513 9.646 
0.501 10.021 
0.463 10.438 
0.482 9.688 
0.485 9.479 
0.507 10.271 
0.546 9.521 


0.529 10.146 
0.532 10.396 


0.474 9.979 
0.499 10.188 
0.490 10.229 
0.451 10.313 
0.496 9.563 
0.468 9.063 
0.524 9.438 
0.526 10.063 
0.521 9.104 
0.465 9.771 


0.457 10.479 
0.535 10.104 


0.488 9.188 
0.515 9.354 
0.543 9.729 
0.510 9.813 
0.454 9.146 
0.479 9.604 


87.582 
85.274 
104.314 
84.632 
96.094 
79.923 
74.388 
70.514 
95.636 
85.478 
68.795 
67.741 
90.679 
77.036 
71.814 
97.981 
70.325 
79.176 
74.470 
71.423 
97.230 
67.129 
85.644 
76.460 
78.493 
70.461 
118.836 


78.170 
73.190 
164.118 
83.386 
84.608 
79.896 
67.622 
67.236 
120.241 
91.709 
55.558 
44.310 
96.088 
84.775 
52.974 
115.708 
72.299 
62.089 
38.660 
66.579 
136.076 
52.081 
113.930 
79.474 
62.743 
44.578 
123.310 


250.555 
127.772 
180.201 
118.151 
125.575 
127.995 
115.862 
113.680 
143.465 
140.391 
174.180 
104.815 
133.996 
116.305 
129.991 
145.906 
117.321 
116.681 
101.617 
120.725 
143.181 
100.040 
144.575 
119.466 
111.810 
98.405 
153.071 


Appendix 2: Experimental design matrix and experimental results 


Table 9. Experimental design matrix and experimental results 


N Design variables Results 

DI H2 R2 D3 R6 D5 D7 E F Stress Stress2 Stress3 
1 13.25 4.86 0.92 12.95 046 1161 11.83 1.93x10? 1101.25 65.68 81.42 389.75 
2 13.50 5.18 1.07 12.68 0.55 11.10 11.59 2.01х10° 1228.75 73.93 71.44 168.58 
3 1296 525 0.91 1219 0.47 11.02 11.91 1.91х105 1206.25 95.50 85.70 181.00 
4 13.38 4.74 1.01 12.88 0.47 11.11 10.99 2.09х10° 1218.75 71.26 67.24 266.04 
> 13.00 4.62 1.07 12.13 045 1127 11.09 2.16х10° 1161.25 94.16 132.11 160.55 
6 13.87 499 1.09 12.18 0.51 1141 11.65 2.08х10° 1178.75 96.32 144.17 160.76 
Т 3.32 5.07 0.92 12.12 0.47 11.82 11.55 1.98х10° 1118.75 154.56 338.55 220.15 
8 13.86 4.77 0.99 12.55 0.52 1190 11.94 1.91х10° 1168.75 117.22 172.13 183.66 
9 3.09 4.51 0.99 12.23 0.45 11.75 11.20 2.14х107 1361.25 167.36 298.76 222.34 
10 13.00 4.94 1.02 12.83 048 11.36 11.28 2.17х10° 1123.75 85.13 85.60 592.49 
11 3.18 4.75 0.90 12.70 0.52 11.85 11.53 1.97х105 1286.25 99.03 116.96 307.71 
12 13.02 4.57 1.00 12.10 0.53 11.16 11.90 1.93х105 1143.75 93.22 118.55 148.49 
13 321 5.00 0.94 12.73 0.51 10.97 11.26 1.92х105 1273.75 70.36 69.97 290.51 
14 3.55 5.33 0.96 12.03 0.50 11.15 11.14 2.02х10° 1343.75 106.86 109.52 161.00 
15 13.14 5.01 0.96 12.05 0.49 11.01 11.80 1.84х10° 1311.25 83.56 95.32 145.15 
16 3.70 4.65 0.95 12.61 049 1121 11.76 2.15х10° 1211.25 101.60 79.63 131.04 
17 13.67 4.97 1.03 1224 0.49 11.85 11.66 1.94х10° 1303.75 171.80 304.14 244.58 
18 3.46 5.11 1.09 12.85 0.51 11.33 1174 1.92х105 1231.25 74.89 76.50 219.69 
19 13.00 4.78 0.97 12.30 0.50 11.03 11.51 2.07х10° 1398.75 104.65 96.00 229.98 


20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 
51 
52 
53 
54 
55 
56 
57 
58 
59 
60 
61 
62 
63 
64 
65 
66 
67 
68 
69 
70 
7 
72 
73 


12.99 
13.20 
3.71 
13.90 
3.05 
3.58 
13.45 
3.65 
3.45 
3.91 
13.19 
3.61 
13.68 
13.29 
13.72 
13.20 
13.10 
13.80 
3.79 
13.42 
3.92 
3.40 
3.30 
3.56 
3.06 
3.64 
13.57 
3.76 
13:53 
13.12 
12.95 
13.05 
13.62 
13.78 
3.59 
13.25 
3.35 
3.69 
13.60 
3.34 
13.90 
2.97 
13.47 
3.44 
13.85 
13.22 
3.80 
13.50 
3.88 
13.40 
3.89 
13.60 
3.17 
3.84 


5.24 
4.75 
4.88 
4.72 
5.29 
5.47 
5.48 
4.82 
5.00 
5.45 
5.03 
4.58 
4.76 
5.02 
5.45 
5.12 
4.55 
4.96 
5.10 
5.13 
5.17 
4.80 
4.59 
4.85 
5.19 
5.37 
4.90 
4.66 
5.35 
4.67 
5.22 
5.34 
5.50 
5.32 
4.79 
5.49 
4.95 
4.91 
5.27 
5.40 
4.81 
5.35 
4.55 
4.87 
4.61 
5.39 
5.04 
4.83 
5.41 
5.38 
5.25 
4.92 
4.70 
5.08 


1.05 
1.06 
0.95 
0.98 
1.01 
1.07 
1.05 
1.04 
1.00 
0.95 
1.08 
1.04 
1.06 
0.97 
1.10 
0.97 
1.08 
1.08 
0.93 
0.93 
0.91 
0.91 
0.90 
0.95 
1.02 
0.92 
0.96 
0.94 
1.04 
0.95 
0.96 
1.08 
1.06 
0.92 
1.04 
1.10 
0.99 
0.91 
0.91 
1.02 
1.00 
1.09 
0.93 
1.03 
1.07 
0.99 
1.01 
1.03 
1.09 
1.05 
1.10 
1.07 
0.96 
0.97 


12.26 
12.92 
12.05 
12.65 
12.40 
12.32 
11.95 
11.96 
12.45 
12.27 
11.98 
12.60 
12.89 
12.85 
12.53 
12.34 
12.59 
12.87 
12.07 
12.20 
12.75 
12.90 
12.40 
12.02 
12.86 
12.84 
12.36 
12.94 
12.56 
12.17 
12.50 
11.99 
12.72 
12.45 
12.65 
12.70 
12.69 
12.47 
12.25 
12.51 
12.64 
12.54 
12.15 
12.66 
12.38 
12.37 
12.22 
12.58 
12.52 
12.50 
12.00 
12.16 
12.55 
12.39 


0.51 
0.50 
0.47 
0.54 
0.49 
0.46 
0.53 
0.53 
0.50 
0.55 
0.48 
0.52 
0.52 
0.50 
0.51 
0.48 
0.51 
0.49 
0.53 
0.52 
0.49 
0.50 
0.49 
0.49 
0.55 
0.47 
0.53 
0.46 
0.45 
0.52 
0.54 
0.50 
0.52 
0.54 
0.46 
0.52 
0.51 
0.46 
0.49 
0.48 
0.51 
0.47 
0.55 
0.52 
0.53 
0.46 
0.47 
0.50 
0.46 
0.53 
0.51 
0.48 
0.53 
0.47 


11.20 
11.44 
1.47 
11.30 
1.40 
1.80 
11.12 
1.60 
1.45 
1.35 
11.26 
1.76 
11.46 
11.93 
11.31 
11.45 
11.23 
11.50 
1.30 
11.70 
1.88 
1.49 
1.08 
1.92 
1.57 
10.99 
11.89 
1.06 
10.95 
11.22 
11.42 
11.94 
11.71 
11.19 
1.55 
11.52 
1.86 
1.24 
11.43 
1.00 
11.32 
1.38 
11.69 
1.04 
11.13 
11.14 
1.66 
11.65 
1.15 
11.45 
1.05 
11.67 
1.29 
1.73 


10.96 
11.81 
11.25 
11.40 
11.24 
11.30 
11.27 
11.20 
11.45 
11.38 
11.37 
11.62 
11.58 
11.79 
11.75 
11.71 
11.36 
11.03 
11.46 
11.75 
11.04 
11.67 
11.25 
11.05 
11.05 
11.44 
11.12 
11.23 
11.13 
11.40 
11.89 
11.56 
11.10 
11.92 
11.45 
11.10 
11.54 
11.69 
11.55 
11.64 
11.57 
11.41 
11.02 
11.88 
11.49 
11.22 
11.52 
11.65 
11.93 
11.86 
11.95 
11.90 
11.45 
11.87 


1.94х105 
2.01х10° 
1.87х105 
1.99х10° 
1.89х10° 
1.80х10° 
1.85х10° 
2.02х10° 
1.99х10° 
2.06х10° 
1.96х105 
2.14х10° 
1.99x10? 
2.13х10° 
1.93х107 
1.88x10? 
1.96х105 
1.97х10° 
2.04х10° 
2.12х10° 
1.80х10° 
1.86х10° 
1.83х105 
2.10х10° 
2.19х10° 
2.08х10° 
1.87x10? 
1.98x10° 
2.03х10° 
1.99x10? 
1.85х105 
1.95x10° 
1.86х105 
2.10х10° 
1.81х105 
2.08х10° 
2.18х10° 
2.06х10° 
2.06х10° 
1.86х10° 
2.12х10° 
2.18х10° 
2.03х10° 
1.91х105 
1.96х105 
2.04х10° 
2.03х10° 
1.89x10? 
2.13х107 
2.09х107 
2.05х10° 
1.83х105 
2.10х10° 
1.79x 10° 


1166.25 
1263.75 
1146.25 
1253.75 
1328.75 
1236.25 
1183.75 
1148.75 
1250 
1196.25 
1181.25 
1116.25 
1251.25 
1318.75 
1238.75 
1348.75 
1386.25 
1111.25 
1158.75 
1291.25 
1326.25 
1113.75 
1186.25 
1388.75 
1366.25 
1346.25 
1391.25 
1256.25 
1133.75 
1241.25 
1298.75 
1358.75 
1353.75 
1203.75 
1243.75 
1176.25 
1306.25 
1128.75 
1296.25 
1156.25 
1221.25 
1316.25 
1233.75 
1363.75 
1271.25 
1283.75 
1258.75 
1331.25 
1266.25 
1191.25 
1163.75 
1308.75 
1226.25 
1383.75 


90.69 
67.36 
107.22 
91.27 
95.85 
123.60 
103.32 
165.33 
88.78 
90.99 
102.47 
102.12 
86.72 
92.66 
78.62 
107.30 
94.41 
80.50 
106.18 
132.42 
124.43 
67.76 
91.21 
320.24 
82.46 
83.03 
139.38 
71.93 
67.20 
99.09 
97.49 
568.32 
102.55 
89.50 
96.02 
79.44 
100.49 
94.18 
108.21 
72.74 
94.29 
93.39 
127.59 
84.98 
96.33 
94.85 
122.59 
106.61 
81.42 
94.42 
87.05 
142.32 
79.47 
129.91 


90.54 
88.74 
189.21 
81.88 
101.76 
148.44 
96.60 
332.65 
94.99 
86.63 
141.68 
135.00 
75.89 
115.60 
78.10 
110.30 
96.77 
66.96 
129.99 
185.09 
100.78 
71.96 
80.74 
2051.87 
110.24 
63.23 
258.13 
60.22 
61.57 
114.70 
101.00 
4072.76 
95.37 
68.21 
90.52 
90.46 
120.63 
72.17 
103.68 
68.06 
77.33 
106.13 
285.42 
75.21 
107.39 
86.82 
197.27 
114.17 
69.63 
87.70 
98.35 
236.27 
87.73 
174.38 


187.89 
485.62 
171.78 
144.70 
240.73 
165.35 
135.32 
215.14 
142.78 
142.73 
157.56 
144.08 
176.46 
331.62 
122.03 
183.24 
304.19 
131.18 
146.99 
330.75 
161.07 
239.71 
152.56 
519.45 
736.30 
185.71 
208.50 
165.41 
131.49 
154.85 
328.55 
847.39 
170.87 
127.77 
150.04 
247.28 
227.94 
130.80 
151.50 
159.67 
134.94 
349.30 
197.43 
197.95 
144.77 
174.50 
189.58 
163.73 
136.90 
152.84 
128.72 
205.51 
223.83 
181.78 


74 13.13 4.80 1.04 12.31 054 1183 11.60 1.88х105 1356.25 144.24 243.30 199.83 
75 13.51 5.14 0.91 12.00 0.52 11.56 11.61 1.98х105 1333.75 148.03 260.73 204.39 
76 3.41 5.44 1.09 12.35 0.49 1.91 11.31 2.04х10° 1301.25 157.11 201.56 194.92 
77 13.36 5.00 0.94 12.62 048 11.84 11.42 2.01х10° 1193.75 104.63 112.33 184.03 
78 3.16 5.46 0.99 12.11 0.54 1.87 11.73 2.17х10° 1131.25 158.52 381.15 240.04 
79 3.39 4.69 1.06 12.48 0.48 1.10 11.18 1.95x10° 1396.25 95.26 87.58 175.54 


80 13.73 5.09 101 12.91 048 11.59 10.95 2.07х10° 1208.75 80.28 76.26 164.02 
81 3.82 4.73 1.03 12.20 0.51 1.62 11.34 211х105 1151.25 117.85 200.77 183.19 
82 3.74 4.56 1.02 12.57 0.53 10.98 11.39 1.84х105 1201.25 73.08 73.65 114.96 
83 3.43 4.90 1.06 12.33 0.51 1.60 11.15 1.81х105 1381.25 127.03 154.60 173.81 
84 13.70 4.68 1.04 1271 0.54 11.09 11.11 1.97х105 1373.75 92.43 65.72 245.16 
85 3.75 5.16 0.93 12.15 0.54 1.20 11.06 1.99х105 1171.25 85.65 95.75 144.25 
86 13.07 530 1.01 12.81 0.51 11.63 11.77 2.17х10° 1276.25 83.49 106.68 547.89 
87 13.54 5.23 0.92 12.42 050 11.05 11.15 2.16х10° 1368.75 90.69 78.46 139.47 


88 13.31 4.64 1.05 12.30 0.45 11.65 11.60 2.15х10° 1293.75 149.23 204.59 184.65 
89 13.01 4.84 1.09 12.10 0.50 11.72 1143 1.90х10° 1138.75 139.41 296.50 197.25 
90 13.85 5.28 0.93 11.97 0.48 11.54 11.21 1.90х105 1223.75 132.20 234.56 185.81 
91 13.30 4.70 1.01 12.08 0.46 11.70 11.30 2.18х10° 1153.75 139.96 331.49 198.89 
92 3.83 5.31 0.95 12.06 0.47 137 11.16 2.11х105 1393.75 25.98 160.64 174.47 
93 13.15 520 0.98 1229 0.50 1148 11.29 1.95х105 1213.75 108.25 108.96 164.78 
94 3.27 5.05 0.96 12.90 0.54 1.79 11.47 2.15x10? 1278.75 00.90 76.56 381.15 
95 3.26 4.54 1.06 1241 0.46 1.00 11.08 2.00х107 1248.75 86.92 81.70 168.28 
96 3.66 471 093 1225 0.52 1.75 10.97 2.14х10° 1313.75 32.79 243.02 202.39 


97 13.11 5.06 1.00 12.75 0.54 1.18 11.70 1.83х10° 1336.25 88.85 83.07 402.50 
98 3.48 5.43 1.03 12.79 0.48 1.25 11.00 2.05х10° 1288.75 81.35 77.71 209.15 
99 3.63 4.60 1.08 12.14 0.53 1.50 11.00 1.94х105 1198.75 29.75 209.66 181.12 
100 13.49 4.65 0.98 12.46 0.54 10.96 11.19 1.89х10° 1173.75 74.67 71.44 130.57 


101 13.81 4.89 0.94 12.63 0.46 1.74 11.78 2.05х10° 1341.25 03.31 123.52 159.00 
102 13.35 4.98 0.98 12.44 0.46 11.90 11.85 1.82х10° 1106.25 110.38 157.44 167.59 
103 13.08 4.50 1.02 12.04 0.50 11.78 11.70 2.09х10° 1246.25 198.78 597.87 278.00 
104 13.94 4.60 1.05 12.60 0.53 11.34 11.72 2.11х105 1376.25 112.22 111.02 152.34 
105 13.10 5.05 0.94 12.43 0.55 11.55 10.98 1.84х10° 1121.25 83.44 97.77 196.20 
106 13.52 5.21 0.97 12.01 0.48 11.35 11.35 1.88х105 1268.75 118.82 149.98 173.65 
107 13.93 5.20 1.03 12.35 0.45 11.53 11.32 2.13х10° 1141.25 92.89 111.00 149.01 


108 12.98 5.42 0.98 12.74 0.54 1.95 11.07 2.00х10° 1136.25 125.69 119.00 525.13 
109 13.24 4.95 0.90 12.28 0.45 11.25 11.80 2.12х10° 1188.75 85.14 89.75 144.63 
110 13.15 4.93 1.07 12.67 0.47 1.68 11.85 2.02х10° 1338.75 100.57 111.93 311.31 
111 3.04 5.40 1.00 12.77 0.49 140 11.17 1.81х10° 1108.75 66.20 83.20 463.98 
112 13.33 4.85 0.97 12.82 0.53 11.58 11.35 1.85х10° 1126.25 74.18 82.46 246.91 
113 13.28 5.10 0.94 12.93 0.47 1.64 11.33 1.82х10° 1323.75 75.48 93.48 406.77 
14 13.37 5.26 0.92 12.80 0.48 11.07 11.63 2.07х10° 1378.75 89.43 71.28 264.36 
115 13.95 5.15 1.02 12.76 0.46 177 11.48 1.80х105 1261.25 97.73 89.91 140.54 
116 13.77 5.15 0.99 12.78 0.50 11.81 11.01 1.87х105 1351.25 95.78 99.06 152.18 
117 13.23 5.30 1.00 12.45 0.54 1.28 11.68 1.90х105 1281.25 79.78 89.59 191.14 
118 13.75 5.36 0.98 12.49 0.55 11.39 11.82 2.16х10° 1216.25 85.45 77.25 134.50 


119 13.45 4.63 1.08 12.21 0.52 11.51 11.50 2.00х10° 1321.25 123.43 200.83 182.56 
120 13.65 4.53 1.05 12.80 0.49 11.80 11.50 1.92x10? 1103.75 92.87 104.34 146.86 
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